home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / DSICO.z / DSICO
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDSSSSIIIICCCCOOOO((((3333FFFF))))                                                            DDDDSSSSIIIICCCCOOOO((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DSICO   - DSICO factors a double precision symmetric matrix by
  10.      elimination with symmetric pivoting and estimates the condition of the
  11.      matrix.
  12.  
  13.      If  RCOND  is not needed, DSIFA is slightly faster.  To solve  A*X = B ,
  14.      follow DSICO by DSISL.  To compute  INVERSE(A)*C , follow DSICO by DSISL.
  15.      To compute  INVERSE(A) , follow DSICO by DSIDI.  To compute
  16.      DETERMINANT(A) , follow DSICO by DSIDI.  To compute  INERTIA(A), follow
  17.      DSICO by DSIDI.
  18.  
  19.  
  20. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  21.       SUBROUTINE DSICO(A,LDA,N,KPVT,RCOND,Z)
  22.  
  23. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  24.      On Entry
  25.  
  26.      AAAA DOUBLE PRECISION(LDA, N)
  27.         the symmetric matrix to be factored.
  28.         Only the diagonal and upper triangle are used.
  29.  
  30.      LLLLDDDDAAAA INTEGER
  31.         the leading dimension of the array  A .
  32.  
  33.      NNNN INTEGER
  34.         the order of the matrix  A .  Output
  35.  
  36.      AAAA a block diagonal matrix and the multipliers which
  37.         were used to obtain it.
  38.         The factorization can be written  A = U*D*TRANS(U)
  39.         where  U  is a product of permutation and unit
  40.         upper triangular matrices, TRANS(U) is the
  41.         transpose of  U , and  D  is block diagonal
  42.         with 1 by 1 and 2 by 2 blocks.
  43.  
  44.      KKKKPPPPVVVVTTTT INTEGER(N)
  45.         an integer vector of pivot indices.
  46.  
  47.      RRRRCCCCOOOONNNNDDDD DOUBLE PRECISION
  48.         an estimate of the reciprocal condition of  A .
  49.         For the system  A*X = B , relative perturbations
  50.         in  A  and  B  of size  EPSILON  may cause
  51.         relative perturbations in  X  of size  EPSILON/RCOND .
  52.         If  RCOND  is so small that the logical expression
  53.         1.0 + RCOND .EQ. 1.0
  54.         is true, then  A  may be singular to working
  55.         precision.  In particular,  RCOND  is zero  if
  56.         exact singularity is detected or the estimate
  57.         underflows.
  58.  
  59.      ZZZZ DOUBLE PRECISION(N)
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDSSSSIIIICCCCOOOO((((3333FFFF))))                                                            DDDDSSSSIIIICCCCOOOO((((3333FFFF))))
  71.  
  72.  
  73.  
  74.         a work vector whose contents are usually unimportant.
  75.         If  A  is close to a singular matrix, then  Z  is
  76.         an approximate null vector in the sense that
  77.         NORM(A*Z) = RCOND*NORM(A)*NORM(Z) .  LINPACK.  This version dated
  78.      08/14/78 .  Cleve Moler, University of New Mexico, Argonne National Lab.
  79.      Subroutines and Functions LINPACK DSIFA BLAS DAXPY,DDOT,DSCAL,DASUM
  80.      Fortran DABS,DMAX1,IABS,DSIGN
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.